The Kondo state in quantum point contacts and the local moment in semiconductor 
quantum dots - two sides of the same phenomenon. 
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This is a three step work: i) we explain why quantum point contacts are similar to ballistic 
quantum dots; ii) we introduce the virtual Kondo state in both systems; iii-lst) this state explains 
0.7 structure in point contacts; iii-2nd) formation of the local moment on this state is described by 
the nearly symmetric Anderson model, we solve it for finite size system having in mind quantum 
dots. We found one large level spacing A* cx (UF) 1 ^ 2 3> A, where U is the charging energy of 
the virtual state, V is the spectral width of this state and A is the mean level spacing of whole 
system. The theory explains periodicity of abnormal level spacing vs gate potential. The theory is 
in agreement with many experiments. 



I. WHY QUANTUM POINT CONTACTS ARE 
SIMILAR TO CHAOTIC QUANTUM DOTS. 

In this work we develop common view on two types of 
experiments, fifSt, differential conductivity of quantum 
point contactsETTJ, and second, fluctuations of level spac- 
ing in semiconductor quantum dotscftr ' . The challenge 
in the former case is to understand the plato of differ- 
ential conductance dl/dV RtJ).7(2e 2 //i), which indicates 
presence of the Kondo statetl. The Kondo state can ap- 
pear from a magnetic impurity, as in the case of metallic 
point contactalj, however the point contact .was made by 
electrodes on top of clean 2d electron gasotJ. 

The level spacing fluctuation in the latter case ia- 
dicates presence of a local moment inside the dotliil. 
There are two possibilities: either magnetic moment is 
formed on the localized Kondo state or the magnetic mo- 
ment isJcpjned due to the Stoner instability of electron 
systemE3~tS ? . However, the clean semiconductor quan- 
tum dots do not have magnetic impurities and formation 
of the Kondo state is unrealistic. The Stoner instability 
occurs undet-Mety restrictive condition met only in some 
experimentaljO. 

The anticipated Kondo state can be found in both sys- 
tems, quantum point contact and semiconductor quan- 
tum dot, as a virtual state bouncing along shortest axis 
of the system. This state is marked by dashed line in Fig. 
fa for contact geometry andiapEig. fb for real potential 
of rectangular quantum dotliaH. We assume that the 
system is ballistic, I £ l p , where / is the mean free path 
and l p is the length of the virtual state. This state is 
very unstable indeed, the electron lifetime on such orbit 
%/Y is of the order of the time of flight l p /vp, in other 
words level uncertainty of such state is of the order of 
level spacing T ~ irTk kp /ml p , 

where kp is the Fermi momentum and m is the ef- 
fective mass. This state is important for two main rea- 
sons i) it-has. good coupling to all other "chaotic" wave 
functionsE2Eil ii) it cost extra energy to put two elec- 
trons on this state compare to all other "chaotic" wave 



functions spread across whole system; the enhancement 
factor is approximately log(kpl p ). 



(a) 



(b) 




FIG. I. Local moment (Kondo state) is formed on unsta- 
ble periodic trajectory between tips of quantum point contact 
(a) or across quantum dot (b). This orbit traps electrons for 
some time and have high charging energy. 

The sweep of the gate potential in both dots and con- 
tacts raise the Fermi level; the -singly occupied Kondo 
state floats on the Fermi levcOEl In this way conduc- 
tance of the point contact at zero bias is given by resonant 
conductance of Kondo state. The level spacing of the 
quantum dot is almost the same as of a non-interacting 
system. When the Fermi level reaches "symmetric" limit 
position, the virtual state can be doubly occupied. The 
system becomes effectively non-interacting. At this gate 
potential one observes kink of the level spacing in the 



1 



quantum dot. In other words, the termodynamic limit of 
the symmetric Anderson model is incompressible electron 
liquid, and finite size system has just pseudogap. 

In the next section we derive effective Hamiltonian by 
making use of Green function technique and arrive at 
the Anderson model. The nearly symmetric limit is bet- 
ter to explore by making use of the Wiegmann solution, 
that is done in Section III. We solve Bethe equation with 
one energy level split out of the band, the "gap" energy 
is exactly the experimentally observed kink in the level 
spacing. Section IV provides analysis of some experimen- 
tal data and summary. 



II. HAMILTONIAN OF THE VIRTUAL STATE. 

Single electron wave functions ipj(r) of a classically 
chaotic system are inhomogeneous. Consider an unstable 
periodic orbit p in a system with two degrees of freedom. 
Such an orbit contributes to the density of states at the 
energy e the sumEa 
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where T p = dS p /de is the period of the orbit, S p (e) is the 
action of the orbit, S p (e) includes the phase changes at 
the conjugated points, A p > is the stability exponent. 

Let us assume that p is the only one relatively stable 
A p ^ 1 periodic orbit in the system. In this particu- 
lar case p contributes a sequence of Lorentzianso to the 
density of states: 
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[S p (e) - 2irnh} 2 + {A p fi/2f 
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Each term in the sum can be regarded as the spectral 
function of the "scar" state 



A n (e) 
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where T n = A p h/(2T p (e n )). 

This spectral function tell us that it is possible 
to construct a non-stationary solution ip n (r) to the 
Schroedinger equation in thepidcinity of orbit p. This 
solution, let us call it "scar"E3, would decay to "flat" 
states with the rate T n /h. The matrix elements of the 
decay process Vj n are not known, but it is clear that 



(4) 



where Ej are energies of the exact wave functions. Fortu- 
nately we will not need explicit values of Vj n in further 
calculations. 



In what follows let us assume that electrons fill the 
system up to the Fermi energy at zero temperature. Let 
n = d is the first "scar" state below the Fermi energy Ep. 
The Hamiltonian for states in the vicinity of the Fermi 
level is 

.Sj<D d+1 

W = El E [ £ i a jv a j° + V oda) a a dlJ + h.< 



<t K j,Sj>D d 

D d = (e d + Cd_i)/2 , 



(5) 
(6) 



and formulation of the Anderson impurity model!! is ac- 
complished by adding the correlation energy 
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with 



n - r 2 



(7) 



(8) 



It should be much larger than the charging energy of the 
dot, Uj, given in general by 



\n - r 2 \ 



(9) 



and almost the same for all flat states ipj. 

The contribution of p to the Green function of the non- 
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where f = (x,y), x is the coordinate along the orbit 
p, the y axis is perpendicular to the orbit p at the point 
x, and W(x) is connected with the second derivatives 
of the action. Then, the width of the "scar" state is 
estimated as ~ yJ%/W. It is clear that W{x) ~ k/l p , 
where k = ^/2med/h, and m is the mass of the particle, 
see for example the analysis of chaotic billiardsEH. The 
integral Eq. (^) for the wave function concentrated in 
the rectangle of the width 



interacting system has the form G(r, r) 



1). 



y/lp/k and length l p gives 
For experiment of Sivan et 
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kl p ~ 100 and the correlation energy U is large. 

To summarize the section both a quantum dot and 
a point contact can have virtual quasi one-dimensional 
state at the energy ed', the charging energy of the state 
is U. This state state is well coupled to a finite num- 
ber of "flat" chaotic states in certain energy interval. So 
we have all ingredients of the Anderson impurity model. 
It allows to compute valence of the impurity state and 
other thermodynamic properties. Our primary interest 
is the compressibility of the electron system at forma- 
tion/disappearance of magnetic moment. This regime is 
described by the almost symmetric limit of the model, 
see the next section. 
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FIG. 2. The function x(A) as it appears in Bethe equations, Eq. |3(j. The rapidities X a lie in the interval [Ao,Am], 
corresponding energies lie between band bottom D and chemical potential fi ~ x(Ao). In the symmetric limit fi ~ + U/2 
and in the assymetric limit n ~ e,j, where e,j is the energy of the virtual state, and U is the charging energy. 



III. THE LEVEL SPACING IN PRESENCE OF 
THE KONDO STATE. 

The purpose of this section is to compute the ground 
state energy En of TV-particle chaotic system in pres- 
ence of the Kondo interaction with quasi-onedimcnsional 
state. Then we will compute the chemical potential /i W = 
En+i — En and the inverse compressibility (the level 
spacing) of whole system A* = /Ujv — /ijv-i- This quan- 
tity is interesting in vicinity of the symmetric regime, 
when spin is just disappears and we expect strong fluc- 
tuations. 

Wiegmann have solved this problem in one 
dimensioned. However dimensionality of the system is 
not important since j-states do not interact with them- 
selves. In order to apply the Wiegmann solution we need 
to postulate the constant energy spacing 2A between j- 
states. This approximation is not bad in chaotic systems 
due to level repulsion. 

At the zero magnetic field the ground state is described 
by the set of rapidities {X a }^ =0 , where M + 1 is number 
of down spins. The rapidities satisfy 

M 



describe the state of the system with N — 2(M + 1) par- 
ticles. The energy of the system is just 



E = 2^x{X a ) 



(12) 



The ground state corresponds to subsequent values of J ai 
fixed by the condition 

J a < Jm , X a < X M , x(X a ) > x(X M ) = D (13) 

where D is the bottom of the energy band, see Fig. |^. 

The ground state is a bit more complicated for odd 
number of particles, N — 2M +1. In this case a = 
Q . . . M — 1, and we introduce new rapidity A: 

1 M ~ 1 1 
J a H — > arctan(A Q — Xg) H — arctan(2A Q — 2A) 

7T ' IT 



-x(X a ) 



5(A) = e d + U/2 - {2UTX) 1 /' 2 = D , 



I AV± -y M-l 

J« + -^arctan(A ct -A /3 ) = --[a;(A Q ) + 0(A)] (10) E = x{X) + 2 ^ x(X a ) . 



(14) 
(15) 

(16) 



13=0 



a=0 



x{X) =e d + U/2 - (C/r) 1 / 2 (A + (A 2 + 1/4) 1 / 2 ) 1 / 2 , (11) 

r = | | 2 is also constant. We keep here terms in the low- 
est order in A. The integer numbers { J a }^ = o completely 



In order to compute the chemical potential we need both 
solutions Eq. ( [l2| ) and Eq. ([l6]), because one additional 
electron changes number of particles from odd to even or 
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vise versa. Fortunately, the difference between the odd- 
electron solution and the even-electron solution is small, 
something happens at the bottom of the band, and not 
near the Fermi surface. For this very reason we will de- 
fine the chemical potential as half of energy one needs 
to add two particles to a system with even number of 
particles: 
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A, = ^f^ 2 . (17) 



Here A» is the level spacing of interacting system. 

When the level spacing of non-interacting system is 
small A — > 0, one can solve Eq. ( |l0| ) by making use of 
the continuous approximation. It is justified by the con- 
dition 



|A +i - A a | < 1 Va . 
So we replace sum Eq. (JlQ) by integral 



(18) 



1 f M iii 1 

J a + — / arctan(A Q - A )cr(A )dX = --rx(X a ) , 

(19) 

where <j(A) ~ A ~~r^ ■ Taking derivative with respect 
to A Q we arrive at the integral equation 



aydX' 



Wa„ 1 + (A-A') 2 



1 dx 
it d\ 



(20) 



which can be formally solved in whole range A e] — 
co, Am] ■ For large Xm we can write equation with a = M 
as Jm + M/2 ~ —D/A, then all J a s become fixed 
J a = -D/A - 3M/2 + a. Integrating Eq. @ from 
— oo to X a (again Am is large) we get equation for X a : 

7T [ " a x dX = e d + U/2-D- 2MA + aA . (21) 

J — OO 

where we made use of x{— oo) = ea + U /2. 
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FIG. 3. The compressibility of electron gas in presence of 
the Anderson impurity. Electron system becomes incompress- 
ible in the mixed- valence regime. 



The general (Wienner-Hopf) solution of Eq. (|2(]) is 
available in the Wiegmann papeiO. The particular case 
Ao = — oo and Am = oo can be solved by Fourier trans- 
form method and the result is 



dk 



1 



2tt cosh(7rA-7rfc 2 /(2C/r)) 



(22) 



it changes a little when Ao,Am are finite but large. For 
large negative A the solution behaves exponentially 



J-oc 2tt 



(23) 



Then we can compute the lowest rapidity Ao from 
Eq. © 



Ao 



1 e d + U/2 - D 
-log — 

7T 



2MA 



V2UT 



(24) 



and indeed it goes to — oo in the symmetric limit of the 
Anderson model 

e d + U/2 — D — 2MA — > 0, ( it is enough to have 
\e d + U/2 + D- 2MA\ < VUf). 

Solution E q. (22| ) works for any Ao if A < Ao- Then we 
combine Eq. (£2) with Eq. ( ^6|) and obtain compressibil- 
ity 
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dN _ 1 

~dfi ~ A J_ QO 2tt cosh(-(7r 2 /4 - jl A )/2jl 2 - k 2 /2) 



(25) 



where fx = (fj, — e d — U/2)/ ^/UT/tt. The compressibility 
is zero at symmetric limit, see Fig. |[ 

In the vicinity of symmetric limit one has simple ex- 
pressions for chemical potential and level spacing 



fi N = ar(Ao) 
= e d + U/2 - 



ttUT/8 



dpi 
dN 



\og(V2UT/A N ) 
t:UT/32 



A N y log 3 (\/2l7r/AAr) 
A N = e d + U/2 - D - (N - 2)A . 



(26) 
(27) 



The level spacing A* diverges for any fixed A in the vicin- 
ity of the symmetric limit. This result indicates that sys- 
tem is unstable near ji = e d + U/2, and experimentally 
one observes kinks of level spacing. 

The kink of the level spacing means that the system 
becomes non-interacting, the chemical potential get its 
U = value 



UN |c/=o = NA + D 
and the observed kink is 



(28) 
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A* — hn\u=o — 



ttUT/8 



V iog(y2c7r/A Ar _ 2 ) 

The maximal possible kink of the level spacing is 



(29) 



' ttUT/4 
log(C/r/2A) 



> A. 



(30) 



It goes to zero in thermodynamic limit A — > 0, but re- 
mains large. We will derive Eqs. (|30|,eq:model.res2) in 
the rest of the section. 

The problem is that the condition Eq. ( |l8| ) is broken 
for the lowest A Q near the symmetric limit. Since for 
Ao — * — oo the density of the rapidi ties goes to zero expo- 
nentially, so |Ai — Ao | oc e 7 ^"' A/y/UT can be very large 
for any fixed A. The main idea of the present paper is 
to introduce one rapidity Ao separated from all the rest 
Ai . . . A a/. We should assume |A - Ai| > 1 in Eq. (0) 
with a = 0: 
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where we expanded arctan(). The integral equation for 
region of " dense" rapidities becomes 
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CTydX' 



1 + (A - A ) 2 tt J Xl 1 + (A - A') 2 
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(32) 



In the lowest order in 1/Ao we neglect the integral in 
Eq. @ and obtain 2c(AqVA = M/2 - J = 2M + D 
that is equivalent to Eq. (pOj). Lowest rapidities are 



Ao m -UT/8A% , 



Ai = 



i log (A N - 2 /y/2U?) , (33) 



where Ai goes to — oo much slowly than Ao. Within the 
approximation made here Ai of system of N + 2 parti- 
cles coincides with Ao of s yst em of N particles. Then the 
obtained pseudogap Eq. ( p0| ) is just x(Xo) — a:(Ai). 

In conclusion of the section we have found the anoma- 
lously large spacing between energy levels of electron gas 
interacting with impurity. This is a finite size effect, 
which requires special treatment. The pseudogap ap- 
pears at the Fermi energy when the model approaches 
complete particle-hole symmetry, /i ss + U /2. 



IV. EXPERIMENTAL DATA AND 
CONCLUSIONS. 

Let us consider experiment, where one fills the system 
by particles. The Fermi level goes up fi(N) ~ JVA + D, 
where N is the number of particles inside the system, 



and D is the band offset. At the value N = Nd given by 
equation 



n(N d ) = e d + U/2 



(34) 



we arrive at the symmetric Anderson model. Just be- 
low this value of N the Fermi energy makes a large jump 
given by Eq. ( |29| ) . Later we will measure the A* in terms 
of mean level spacing 



<5* = A*/A « y/(n/4)x/\ogx , x = UT/2A 2 > 1 . 

(35) 

The Fermi energy will jump again near N = Nd+i- The 
distance between jumps is Nd+i — Nd ~ {e d +i — ed)/ A is 
inverse proportional to the length of the orbit l p 



2L 



SN = N d+1 -N d ~ — v/27rA d+ 



(36) 



This formula was derived for A = ttH /(mA), where 
L is size of the quantum dot, A m L 2 is its area and 
SN < N d+1 . 

The experimental data from few references is summa- 
rized in Table |. The first block contains data taken from 
cited papers, then we provide estimation of kp, A, T, U, 
SN, and <5*. Last block in Table | contains SN and <5» 
taken from charts of -E/v+i — En vs N reported in cited 
papers. 

Typically l p ~ 2L and we will use SN — y/2nNd+i. 
The Coulomb energy of the virtual state is U ~ 
Ec^ogkplp where Ec is the charging energy of the dot, 
and kp is wave number on Fermi surface. The coupling 
r is of the order of level spacing of virtual state, at least 
it larger than level spacing of host system: A < T < 
e d +\ — ed- For this reason Table || contains two values 
for r - minimal and maximal estimates. Last computed 
value is 5*, which eventually have minimal and maximal 
estimates derived from Eq. (BflLAVe see that our theory 
works well for clean samplesoofi. and it is not applica- 
ble for disordered quantum dotsa. The virtual state is 
absent in perfectly rectangular quantum dotEl (without 
shape deformations) as well. 

The conductance of a quantum point contact (QPC) 
with the. Kondo state inside has been computed 
recentljcil and the model explained all the experimen- 
tal data. Existence of the Kondo state in a QPC is very 
plausible, and therefore we expect to see Kondo state in- 
side quantum dots too. The present theory should not 
be confused with Kondo effect due to charging of entire 
doiH. However the intersting feature of both types of 
experiments (Kondo state on "scar" and Kondo state in- 
side entire dot) is that the sate voltage allows to scan all 
regimes of Anderson mode£j. 

To summarize, inhomogeneities of chaotic wave func- 
tions can be gathered together into an additional state 
of small size. In this way one arrives at Anderson's im- 
purity model with finite number of free electron states. 
At certain value of the chemical potential the model has 
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complete particle-hole symmetry. The chemical potential 
does not approach this value gradually, but rather irregu- 
larly. The energy scale of this effect is a new combination 
of the parameters of Anderson's model and the mean level 
spacing. This scenario is a possible explanation of the ir- 
regularities observed experimentally in quantum dots. 
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TABLE I. Measured data / Intermediate calculations/ 
Observed values of SN and <5». energies are given in meV, 
length are given in fim and wavenumbers in /urn -1 . 
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0.02 
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0.05 - 2.0 
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0.02 - 1.0 
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25 


46 


5* 


3 - 10 


4 - 16 


3 - 7 


3 - 13 
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15 
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<5* 


10 


6-8 
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f - circumference J - the sample # 1 
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